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Abstract: Three-particle correlations in quark and gluon jets are computed for the first time in perturba- 
tive QCD. We give results in the double logarithmic approximation and the modified leading logarithmic 
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bution. We thus provide a further test of the local parton hadron duality and make predictions for the 
LHC. 
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1 Introduction 



The observation of quark and gluon jets has played a crucial role in establishing QCD as the theory of 
strong interaction within the standard model of particle physics. The jets, narrowly coUimated bundles of 
hadrons, reflect configurations of quarks and gluons at short distances. Powerful schemes, like the double 
logarithmic approximation (DLA) and the modified leading logarithmic approximation (MLLA), which 
allow for the perturbative resummation of soft-collinear and hard-coUinear gluons before the hadroniza- 
tion occurs, have been developed over the past 30 years (for a review see lH). One of the most striking 
predictions of perturbative QCD, which follows as a consequence of angular ordering (AO) within the 
MLLA and the local parton hadron duality (LPHD) hypothesis [2], is the existence of the hump-backed 
shape H] of the inclusive energy distribution of hadrons, later confirmed by experiments at colliders. 
Indeed, the shape and normalization of single inclusive distributions are compared with an experiment; 
a constant /C^'^, which normalizes the number of soft gluons to the number of charged detected hadrons 
(mostly pions and kaons), turns out to be close to unity (JC^^ ~ 1), giving support to the similarity be- 
tween parton and hadron spectra ||T]. Thus, the study of inclusive observables like the inclusive energy 
distribution and the transverse momentum k± spectra of hadrons |I3] has shown that the perturbative stage 
of the process, which evolves from the hard scale or leading parton virtuality Q ~ £^ to the hadronization 
scale Qq, is dominant. In particular, these issues suggest that the hadronization stage of the QCD cascade 
plays a subleading role and, therefore, that the LPHD hypothesis is successful while treating one-particle 
inclusive observables. 

The study of particle correlations in intrajet cascades, which are less inclusive observables, focuses 
on providing tests of the partonic dynamics and the LPHD. In [4J, this observable was computed for 
the first time at small x (energy fraction of the jet virtuality taken away by one parton) in MLLA for 
particles staying close to the maximum of the one-particle distribution. In [5J, the previous solutions 
were extended, at MLLA, to all possible values of x by exactly solving the QCD evolution equations. 
This observable was measured by the OPAL Collaboration in e^e^ annihilation at the peak, that is, 
for ^/s = 91.2 GeV at LEP [6|. Though the agreement with predictions presented in 15] was improved, 
a discrepancy still subsists pointing out a possible failure of the LPHD for less inclusive observables. 
However, these measurements were redone by the CDF Collaboration in pp collisions at the Tevatron for 
mixed samples of quark and gluon jets Q. The agreement with predictions presented in ||4l turned out 
to be rather good, especially for particles having very close energy fractions (xi X2). A discrepancy 
between the OPAL and CDF analysis showed up and still stays unclear. Therefore, the measurement of 
the two-particle correlations at higher energies at the LHC becomes crucial. Furthermore, going one step 
beyond, in this article we give predictions for the three-particle correlations inside quark and gluon jets. 
This observable and the two-particle correlations can be measured on equal footing at the LHC so as to 
provide further verifications of the LPHD for less inclusive observables. 

2 Kinematics and evolution equations 

A generating functional Z{E, 0; {u}) can be constructed |[T] that describes the azimuth averaged parton 
content of a jet of energy E with a given opening half-angle ©; by virtue of the exact AO (MLLA), 
which satisfies an integro-differential system of evolution equations. In order to obtain exclusive n- 
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particle distributions D^^\ki, E) one takes n variational derivatives of Za over u{ki) with appropriate 
particle momenta, i = 1 . . . n, and sets u = afterwards; inclusive distributions are generated by 
taking variational derivatives around u = I. Let us introduce the n-particle differential correlations for 
A = G, Q, Q jets as, 

^S"U^) = - • • • -D^A^ f a . . . ^, In (1) 

Z Z \ Z Z Qq J 

together with = ^\].rSX) for later use; Xi corresponds to the Feynman energy fraction of the jet 

taken away by one particle "i" and z is the energy fraction of the intermediate parton. For instance, for 
three-particle correlations n = 3, the observable to be measured reads C)^^^^ = ^^^^^3 ■ The production 
of three hadrons is displayed in Fig|T] after a quark or a gluon {A) jet of energy E with half opening 
angle 0o and virtuality Q = EQq has been produced in a high energy collision. The kinematical 
variable characterizing the process is given by the transverse momentum k_i = zEQi > Qo [or (1 — 
z)EQi > Qo] of the first splitting A — BC. The parton C fragments into three offspring such that 
three hadrons of energy fractions xi, X2, and X3 can be triggered from the same cascade following the 
condition ©q > 0i > 02 > ©3, which arises from exact AO in MLLA |[1]. We make use of variables. 



0) 

and A = ln{QQ/ Aqcd)- The two variables entering the evolution equations are z and ©i, such that 
xi < z < I ^ < £ < ii. Accordingly, the anomalous dimension related to the coupling constant can 
be parametrized as 



where /3o = 4^ (t-^^ — f'^/^ij)' with Tr = 1/2 and Uf the number of light quark flavors. From 
AO and the initial condition at threshold x^EQq > x^EQi > x^EQ^ > Qq, one has the bounds 
< ©1 < ©0 =^ < y < 7/3 for the integrated evolution equations. The evolution equations satisfied 
by ^ are derived from the MLLA master equation for the generating functional Za{E, 0; u{ki)). For 
three-particle correlations, one takes the first second ^^^(0^, and finally third Su{kf)l5u{k,) 

functional derivatives of Za{E, ©; u{ki)) over the probing functions u{ki) so as to obtain the differential 
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system of evolution equations: 



G?J = llG^'^-^ll {of -/3o7o^gW) + (a - h)^l [(df^G.+Gf^G, + G^^)-^G,)^ (3) 



where A 



(2) 



(2) 



^12 ^3 



ii?^2 



" (2) 

A2'^Ai. The subscripts 



and = A^^) - A,A,A, - r.,, ^3 - -13 

i and y in Eqs. Q and Q denote and d/dy, respectively. The first terms of Eqs. and ^ 
are of classical origin and, therefore, universal. Corrections cx — |, a, (a — 6), and (a — c), which are 
0{y/a~s) suppressed, better account for energy conservation at each vertex of the splitting process, as 
compared with the DLA 0(1). The hard constants are obtained after integration over the regular part of 
the DGLAP splitting functions IB as performed in ||4]|5]. In the equation for the gluon initiated jet Q, 
the first and second constants a{nf = 3) = 0.935 and b{nf = 3) = 0.915 were obtained in the frame of 
the single inclusive distribution and two-particle correlations respectively ||4]. The third constant c(nj) 
appearing for the first time in this frame reads 



c{nf) 
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0.917. 



2.1 MLLA and DLA solutions of the evolution equations 

Equation Q is self-contained and can be solved iteratively by setting G^^^ = Cq^^^^GiG2Gs and 
Gj-j^ = Cq^GiGj in the left- and right-hand sides of Accordingly, the solution of ^ is also 
obtained by setting ^ (jg) Q^Q^Q^ and Q^f = C^a^.QiQj in the left-hand side of © and 
G'(^) = c'^Q^^^G\G2Gz in the right-hand side of the same equation such that the iterative solutions can 
be written in the compact form 

The MLLA two-particle correlators C\^^ will be taken from 15] for the computation of ■ Moreover, 

= 1 + j-^V^^-^^ , (5) 
2 + Ai2 + Ai3 + A23 + ei 



,(3) _ l -C^. + gP+ef -6 1 

2 + A12 + Ai3 + A23 + ei 



G123 o , A , A , A , _ ^ ^ 



and for the quark jet 



^(2) ^ Ci - h 

3 + A12 + Ai3 + A23 - a^e + h ' 

3 + Ai2 + Ai3 + A23-a^'£ + ei ' 

where ^'^ = V'l,^ + "02/ + '>p3/ = 0(70) and -0 = ln[G(£, y)]. Higher order corrections arising from the 
solution of the system of Eqs. [2] and |3] have been neglected in In this case, G{i, y) is the inclusive 



3 



energy distribution, which will be inserted from the steepest descent method presented in [5 |. The other 
functions appearing in ^ and ^ are Aij = 7q"^ (V'i/V'j,j/ + i'i,y'4'j/) = C'(l) and 



C, 
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eg) 

(-'123 
(j(2) 



O(7o), Cy 



c, 



(2) 



o); Xy 



M3) 
Gi23,y 

eg) 

^(2) 



(2) 



lo [X7'^y + Xl''i'e)=0{jo), 



0), 



InC^^) and v 

(-'123 ^ 



with C 
previous by replacing 

^(3) 



in . The set of functions appearing in dTJl and dUl is obtained from the 



eg) andC, 



(3) 



eg) where the dotted C^"^ and 



(2) 



Cj^ are the DLA solutions of the two- and three-particle correlators; that is why this solution is said to 
be iterative. Moreover, corrections ei, ei and Ci ,Ci ^re very small and do not play a significant role in 
the shape and normalization of the three-particle correlations. 

The DLA two-particle correlators are taken from ||8] and the DLA expression for cg), can be obtained 
from ^ by setting all MLLA 0(70) corrections to zero: 



C 



^(3) 



(2) 



1 



CaI + \ 



Nr. 



(2) 



C 



1 

(2) 



^(2) 



^(2) 



c 



A2: 



(9) 
(10) 



I + c 



^(2) 



Ca 



2 + A12 + A13 + A 



+ 



^23 



Ci 2 + A12 + Ai3 + A23 



The solutions have the following simple physical interpretation: the first term (= —1) in the left-hand 
side translates the independent or decorrelated emission of three hadrons in the shower. After inserting 
the two-particle correlator with color factor oc ^ Q in the left-hand side of (ITOl ). terms oc ^ corre- 
spond to the case where two partons are correlated inside the same subjet, while the other one is emitted 
independently from the rest. Next, replacing (|9]l in the right-hand side of (ITOl ). one obtains a contribution 



oc 



„2 describing the independent emission of two partons inside the same subjet. The last term oc „2 
involves three particles strongly correlated inside the same partonic shower as depicted in Fig|T] This 
term is indeed the cumulants of genuine correlations, first obtained in this article for this observable. 

The evaluation of (O, which is expressed in terms of the logarithmic derivatives of the single inclusive 
distribution \d[G{1, y)], will be performed using the steepest descent method to determine y) |l5l[8]|. 
Thus, the MLLA logarithmic derivatives were written in ||5] in the form: 



i^i,yif^,'^) 



A/if^i^ ^i) = 7oe'^' + 77^70 e^^'Qif^i, n) - tanhfj - tanhz^j coth/ijl 1 + e^^'Qim, Ui 



(11) 



7oe 
1 



7^70 



2 + e J^i) + tanh — tanh i/j coth ( 1 + e ^^Qdj-ijV'i) 



/3o7o I + ta.nhud 1 + K{fii,Ui)) -C{fii,Ui){l + e-'''Q{fii,Ui)) +O(7o), (12) 
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where the functions Q{^i,i'i), C{fii, Vi) and K{^i, ui) are defined in |5j and (/ij, fj) are expressed as 
functions of the original variables {£{, yj) by inverting the nonlinear system of equations O: 



Vi- £i _ (sinh - - (sinh 2ui - 2vi) 



+ Vi 2(sinh^ /ij - sinh^ fj) 

sinh i/i sinh fii 



In particular, this method allows for the estimation of the observable for particles with energies near 
the maximum or hump {£max = ^/2) of the one-particle distribution \ i — Y/2 a oc Y^^'^, which 
applied to the three-particle correlations will appear in a forthcoming paper. For instance, at DLA one 
has Ajj = 2cosh(/ij — fXj) with such a parametrization of the logarithmic derivatives of the inclusive 
spectrum. Close to the hump one has Ajj ~ {£i — Ij)'^; thus the correlations are expected to be quadratic 
as a function of {li — £j ) and to have a maximum for particles with the same energy Xi = xj . In this frame, 
the role of MLLA corrections should be expected to be larger than for the two-particle correlations. 
Indeed, higher order corrections increase with the rank of the correlator, which is known from the Koba- 
Nielsen-Olesen problem for intrajet multiplicity fluctuations [91. For the two-particle correlations, for 
instance, one has oc —b{ipi^i + ip2,e) and for the three-particle correlator one has the larger correction 
(X -c(V'i/ + ^2/ + ip3,e)- 

2.2 Phenomenology and comparison with existing e+e~ and pp data 

The study of n-particle correlations is very important because, being defined as the n-particle cross 
section normalized by the product of the single inclusive distribution of each parton 



the resulting observable becomes independent of the constant JC , thus providing a refined test of QCD 
dynamics at the parton level. Since our study of three-particle correlations depends on previous results 
for two-particle correlations, we briefly review recent results about this observable. The MLLA evolution 
equations for two-particle correlations, quite similar to those leading to the hump-backed plateau, were 
solved iteratively in terms of the logarithmic derivatives of G{i, y) [5j. That is how, the result previously 
obtained by Fong and Webber in ||4J, only valid in the vicinity of the maximum £rnax of the distribution, 
was extended to all possible values of x. Consequently, as displayed in Fig|2j the normalization of the 
more accurate solution of the evolution equations is lower and reproduces some features of the OPAL 
data at the peak Q = 91.2 GeV of the e^e~ annihilation, like the flattening of the slopes towards 
smaller values of x JH. Qualitatively, our MLLA expectations agree better with available OPAL data than 
the Fong-Webber predictions ||2]. There remains however a significant discrepancy, markedly at very 
small X. In this region nonperturbative effects are likely to be more pronounced. They may undermine 
the applicability to particle correlations of the LPHD considerations that were successful in translating 
parton level predictions to hadronic observations in the case of more inclusive single particle energy 
spectra {\\. 

These measurements were redone by the CDF Collaboration for pp collisions at center of mass energy 
^/s = 1.96 TeV for mixed samples of quark and gluon jets [7 |. For comparison with CDF data, the 
two-particle correlator was normalized by the corresponding multiplicity correlator of the second rank. 
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Figure 2: Two-particle correlations in two quark jets yR = ^ + j 121 in the process e+e — ^ qq as 
a function of ii + £2 = \ ln(xiX2)| for £1 — £2 = ln(x2/xi) = 1.0 (left) and £1 — £2 = \ ln(x2/2;i)| for 
£1 +^2 = ln(xi3;2) = 6.0 (right). 



which defines the dispersion of the mean average multiplicity inside the jet. In this case, the MLLA 
solution by Fong and Webber ||4l, the more accurate MLLA solution ||51, and the NMLLA solution |[3l 
were compared with the CDF data. The Fong-Webber predictions turned out to be in good agreement 
with CDF data in a range from large to small x, also covering the region of the phase space where MLLA 
predictions should normally not be reliable, that is, for x > 0.1 (see FigjSll. As these figures were taken 
from llll, different notations have been used in this case, for instance, i = ^ = ln(l/x), = ^ — £,max 
i^max = imax = \ HQ/Qo)) such that A^i + A^2 = ^1 + ^2 " HQ/Qo) and A^i - A^2 = h- 
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Figure 3: Two-particle correlations in a mixed sample of gluon and quark jets in collisions as a 
function of A^i + A^2 = | ln(3;i3;2)| - ln(g/Qo) for A^i = A^2 (left) and A^i - A^2 = | ln(x2/xi)| 
for A^i = -A^2 (right). 



As observed in Fig|3](left), the data are well described by the three cases in the interval A,^i + A^2 > 
—0.5, that is, at very small x. However, the Fong and Webber's solution also describes the data for 
A^i + A^2 < —0.5, that is, for larger values of x where the MLLA is no longer valid. QCD color 
coherence for Figj3](left, the peak at A^i + A^2 = —1-5 is due to numerical uncertainties) should be 
observed if the analysis is extended to A^i + A^2 > 2.5. Moreover, the NMLLA solution 131 extends. 



6 



like for the k± spectra, the region of apphcability of such predictions for larger values of x. In Q, it was 
concluded that despite the disagreement with the OPAL data in Fig|2l the LPHD stays successful for the 
description of less inclusive energy-momentum correlations. Therefore, in this paper we encourage the 
analysis of these observables by other collaborations like ALICE, ATLAS, and CMS in order to clarify 
this mismatch. 

Y=7.5 ;^=0.1 ; |Ln(x^X2)|=10 Y=7.5 ;>.=0.1 ; ILnCx^x^jhlO 




Figure 4: Gluon jet 3-particle correlator as a function of | ln(xiX2)| for xi = X2 and ln(l/x3) (left) and 
as a function of ln(a;2/a:;i) for fixed | ln(xia;2)| and ln(l/x3) (right). 



3 Predictions for three-particle correlations and phenomenology 

Finally, in order to extend the applicability of the LPHD to a larger domain of observables, we per- 
form theoretical predictions for three-particle correlations in the limiting spectrum approximation (Qq « 
Aqcd)- This observable and two-particle correlations can be measured on equal footing at the LHC. 
We display the MLLA solutions (01) of the evolution equations ^ and ([3]). The correlators are functions 
of the variables £i, yi and the virtuality of the jet Q = EQq. After setting yi = Y — ii with fixed 
Y = ln((5/Qo) in the arguments of the solutions dUl, the dependence can be reduced to the following: 
cSljh,i2,i3, Y) and C^^ljh,i2,i3, Y). 

In Fig. H the DLA ([lOll and MLLA dH) three-particle correlators for A = G and A = Q , Q, 

^(3) ^ '^123 A3) ^ Vl23 

^^23 G1G2G3' '^'^^ Q1Q2Q3 

are displayed, respectively, as a function of the difference (£1 — £2) = ln(x2/a;i) for two fixed values 
of £3 = ln(l/3;3) = 4.5, 5.5, fixed sum {h + £2) = \ ln(xiX2)| = 10, and, finally, fixed Y = 7.5 
(virtuality Q = 450 GeV and Aqcd = 250 MeV), which is realistic for LHC phenomenology ||5l. The 
representative values £3 = ln(l/x3) = 4.5, 5.5 (X3 = 0.011, X3 = 0.004) have been chosen according 
to the range of the energy fraction Xj <^ 0.1, where the MLLA scheme can only be applied. 

In Fig. m the DLA and MLLA HI) three-particle correlators for A = G and A = Q,Q are depicted, 
in this case as a function of the sum {£1 + £2) = \ ln(xiX2)| for the same values of £3 = ln(l/x3) = 
4.5, 5.5, for xi = X2 and Y = 7.5. As expected in both cases, the DLA and MLLA three-particle 
correlators are larger inside a quark than in a gluon jet. Of course, these plots will be the same and the 
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Figure 5: Quark jet 3-particle correlator as a function of | ln(xiX2)| for xi = X2 and ln(l/x3) (left) and 
as a function of \n{x2/xi) for fixed | ln(xiX2)| and ln(l/x3) (right). 

interpretation will apply to all possible permutations of three particles (123). As remarked above, the 
difference between the DLA and MLLA results is quite important in pointing out that overall corrections 
in 0{^/a^) are large. Indeed, the last behavior is not surprising as it was already observed in the treatment 
of multiphcity fluctuations of the third kind given by = 4.52 [l - (2.280 - 0.018n/)^/al] ifTOl . 

For instance, for one quark jet produced at the peak of the e+e^ annihilation (Q = 45.6 GeV), one has 
as = 0.134. Replacing this value into the previous formula for the quark jet multiplicity correlator, one 
obtains a variation from 4.52 (DLA) to 0.83 (MLLA). Because of this, DLA has been known to provide 
unreliable predictions which should not be compared with experiments. From FigHJ the correlation is 
observed to be the strongest when particles have the same energy and to decrease when one parton is 
harder than the others. Indeed, in this region of the phase space two competing constraints should be 
satisfied: as a consequence of gluon coherence and AO, gluon emission angles should decrease and on 
the other hand, the convergence of the perturbative series k± = XiEQi > Qq should be guaranteed. That 
is why, as the collinear cutoff parameter Qq is reached, gluons are emitted at larger angles and destructive 
interferences with previous emissions occur. This effect is clearly observed in Fig. 51 the steep fall of the 
distribution is more pronounced in the quark jet than in the gluon jet. Moreover, the observable increases 
for softer partons with X3 decreasing, which is for partons less sensitive to the energy balance. In Figl5] 
the MLLA correlations increase for softer partons, then flatten and decrease as a consequence of soft 
gluon coherence, reproducing for three-particle correlations the hump-backed shape of the one-particle 
distribution. Because of the limitation of phase space, one has C^^^ < 1 for harder partons. 

4 Summary 

In this paper we provide the first full perturbative QCD treatment of three-particle correlations in parton 
showers, provide a further test of the LPHD within the limiting spectrum approximation, and briefly 
revise the comparison of two-particle correlations with OPAL and CDF data. The correlations have been 
shown to be strongest for the softest hadrons having the same energy xi = X2 = x-^ in both quark and 
gluon jets, increasing as a function of ln{xi/xj) and | ln(xjXj)| when Xk softens, that is for partons less 
sensitive to the energy balance. This result becomes therefore universal for n-particle correlations. 
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Coherence effects appear when one or two of the partons involved in the process are harder than the 
others, thus reproducing for this observable the hump-backed shape of the one-particle distribution. Also, 
the two- and three-particle correlations vanish (C^^^ — >• 1) when one of the partons becomes very soft, 
thus describing the hump-backed shape of the one-particle distribution. The reason for that is dynamical 
rather than kinematical: radiation of a soft gluon occurs at large angles, which makes the radiation 
coherent and thus insensitive to the internal parton structure of the jet ensemble. 

We give the first analytical predictions of this observable in view of forthcoming measurements by AT- 
LAS, CMS, and ALICE at the LHC. Further information from the comparison with forthcoming data 
may also help to improve Monte Carlo event generators in the soft region of the phase space in intrajet 
cascades, where PYTHIA, ARIADNE and HERWIG face difficulties while reproducing the data ifTll . 
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